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Abstract. We consider a one dimensional transport model with nonlocal velocity given by the 
Hilbert transform and develop a global well-posedness theory of probability measure solutions. Both 
the viscous and non- viscous cases are analyzed. Both in original and in self-similar variables, we ex- 
press the corresponding equations as gradient flows with respect to a free energy functional including 
a singular logarithmic interaction potential. Existence, uniqueness, self-similar asymptotic behavior 
and inviscid limit of solutions are obtained in the space 7^2 (R) of probability measures with finite 
second moments, without any smallness condition. Our results are based on the abstract gradient 
flow theory developed in 2:. An important byproduct of our results is that there is a unique, up 
to invariance and translations, global in time self-similar solution with initial data in 7^2 (R), which 
was already obtained in 17, 6 by different methods. Moreover, this self-similar solution attracts all 
the dynamics in self-similar variables. The crucial monotonicity property of the transport between 
measures in one dimension allows to show that the singular logarithmic potential energy is displace- 
ment convex. We also extend the results to gradient flow equations with negative power-law locally 
integrable interaction potentials. 
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1. Introduction 

In this work, we are interested in developing a well-posedness theory of measure solutions to the 
equation 



n n / ut + {H{u)u)^ = 

^ ' 1 n(x,0)=no(x) 



H{u) = -P.V. [ 
71" Jr 



■ dz . 



with general nonnegative initial Borel measures uq. Here, the term H(u) denotes the classical Hilbert 
transform 

u{z) 

X 

Since the equation is of transport nature and in divergence form, we expect sign preservation and 
mass conservation. Therefore, we will restrict our attention to probability measures as initial data. 
This equation is nothing else than a ID-dimensional continuity equation in which the velocity field 
is given by the Hilbert transform and it has been proposed as a simplified model in fluid mechanics 
[13j and in dislocation dynamics [6] as we will discuss in next subsection. This equation can be 
formally considered a particular example of the theory of gradient flows in the space of probability 
measures [2] as it will be further elaborated in subsection 1.2. 

The main aim of this work is to show that unique measure solutions of gradient-flow type can be 
constructed for the problem 



(1.2) 



u{x, 0) = uo{x) 
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with K > and uq € 7^2 (1^) the set of probabihty measures on the real hne with finite second 
moments. Moreover, the solutions will continuously depend on both the initial data uq and the 
viscosity parameter k > 0. The main tools of this construction are the variational schemes based 
on optimal transportation theory originated for the seminal work [24J. 

Moreover, we will be able to characterize the large time behavior of the solutions. In fact, we 
show that suitable scaled equations related to (jl.ip and ()1.2p have unique stationary solutions fixed 
by the normalization of the mass. Furthermore, we are able to show that the solutions constructed 
converge for large times to these stationary solutions exponentially fast in some transport distance. 
These stationary solutions correspond to self-similar solutions for the original equations. 

This manuscript is organized as follows. In the next two subsections we recall the main results 
already obtained in the literature and the origin of these models. On the other hand, we introduce 
some basic notations and definitions about optimal mass transportation theory essential to our 
construction. Section 2 is devoted to introduce self-similar variables and rewrite our problem in 
the form of a gradient fiow in the space of probability measures of a free energy functional. Key 
properties of the functionals involved are shown in subsections 2.1 and 2.2. Finally, we state and 
prove our existence, asymptotic behavior and inviscid limit results in Section 3. 

1.1. Motivation: Fluid and Fracture Mechanics. One of the motivations to analyze these 
equations arose from the mathematical fluid mechanics literature. In fact, it appears as a simplified 
one dimensional model [3] mimicking the structure of the 3L'-Navier-Stokes equations and the 2D- 
quasi-geostrophic equations p^: 

r ut + {e-V)u = 
(1.3) <^ 9 = S/^(j), u= -(-A)3(/, 

[ u{x,0) = uo{x), 

where V"*" = {—d2,di) and u{x,t) represents the air temperature. Since di'v{R^u) = 0, rewriting 
the system (|1.3p in terms of the Riesz transform given by 

R,{u){x,t) = ^P.V. f ^-^l^u{y,t)dy 
27r J^2 \x-y\-^ 

the result is 

/-I / ut + div [{R^u) u] =0 

^ ' I n(x,0)=no(x), 

Considering u{x, t) with (x, t) € [0, oo) x R and replacing the Riesz transform in (|1.4p by the Hilbert 
transform in one dimension leads to (jl.ip or ()1.2p with diffusion, see [13] for previous related works 
and simplifled models. Mathematical fluid mechanics arguments have been used to analyze existence 
and uniqueness, finite time blow-up of smooth solutions, and other issues, see [27 } [T3 l [T5l [TBI [T8 | 
and the references therein related to these equations and other nonconservative variants. 

More precisely in our case, sign-changing periodic C^-solutions of (|l.ip blow up in finite time, 
in the sense that its C^-norm diverges in finite time as shown in [13j. On the other hand, global 
existence and uniqueness of smooth solutions for the Cauchy problem on the whole real line is proved 
in [12] for strictly positive initial data for (jl.ip and for general nonnegative initial data for (|1.2p . 
The same authors show that for non negative touching-down initial data the Cauchy problem for 
p.ip is locally well-posed for smooth solutions and that solutions do blow up in finite time in the 
norm. 

Apart from the structural similarities, the equations (|1.4p and (jl.ip have different properties. For 
instance, while the first one has a Hamiltonian structure, the second one being one dimensional can 
be considered rather as a gradient fiow as we will discuss in next subsection. 

The other source of motivation to analyze equations (jl.ip and (jl.2p comes from dislocation dy- 
namics in crystals |20i I21j [22 l I17 j . Here, the unknown u represents the number density of fractures 
per unit length in the material. The existence of explicit self-similar solutions and the convergence 
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towards them was studied in [n\ |6] showing that nonnegative solutions play an important role in 
the large time asymptotics of p.2p and related problems. In fact, we will give a characterization of 
the self-similar solution as the minimizer of a free energy functional intimately related to its gradient 
flow structure. In this way, we will show that the solution does really converge in suitable scaling 
and in transport distances to the self-similar profile. 

1.2. Gradient Flows for Probability Measures. Let us remind some basic facts about optimal 
mass transport, which will be useful to our study of solutions of the Cauchy problem ()1.2p . For 
more details we refer the reader to [301 [2]. Let us denote by 'P(M'^) the space of probability measures 
on W^. We start reminding the definition of push forward of a measure p G V{M.'^). 

Definition 1.1. Let p be a probability measure on R'^ and let T : — >■ M 6e a Borel map. The push 
forward T^p G V{W^) of p through T is defined by T^p{I) := p(T~^{I)) for any Borel subset / C M. 
The measures p and T^p satisfies 



f{T{x))dp{x)= / f{y)dT^p{y), 

for every bounded or positive continuous function f . 

Let /O, /i S 'P(M'^) and Tp : M'^ — )• M"^ such that T^^p = p. The map Tp is called a transport map 
between the probability measure p and p. We also recall the notion of transport plan between two 
probability measures. 



Definition 1.2. Given two measures p and p ofV{M. ) the set of transport plans between them is 
defined by 

r(p, p) := {7 e X m ■■ ^« 7 = p, = A ' 

where vr* : M'^ x M*^ M'^, i = 1, 2 are the projections onto the first and second coordinate: 7^^{x, y) = 
X, Tr'^{x,y) = y. In other words, transport plans are those having marginals p and p. 

Our aim is to study solutions of the Cauchy problem (jl.ip and ()1.2p in an appropriate subspace 
of 7-'(M) endowed with a transport distance, the so-called euclidean Wasserstein distance. Consider 
the set 

p2(M'^) = |p € V{M.'^) : j \x\'^dp{x) < ooj . 
The euclidean Wasserstein distance is defined on 7^2 (I^'^) as: 

Definition 1.3. For any probability measure p, p £ 7^2(1^'^) the euclidean Wasserstein distance 
between them is defined by 




d2{p, p) := min { { / \y - x\ dj{x,y) ] ■.-feT{p,p) 



We denote by Tq{p,p) the set of optimal plans, i.e., the subset of T{p,p) where the minimum is 
attained, i.e, 

ro(p, = i 7 e T{p, p): I \y- x\'^d'y{x, y) = dl{p, p) \ . 

The space 7^2 (1^'^) endowed with d2 becomes a complete metric space. The convergence in d2 is 
equivalent to weak-* convergence as measures together with convergence of the second moments, 
see |3m Theorem 7.12]. We will denote by 'P2'^{W^) the subset of probability measures with absolutely 
continuous densities with respect to Lebesgue measure and finite second moments. It is well-known 
that for any p, p ^ V2^{M.'^), the minimum in the definition of d2 is achieved by a plan defined by 
an optimal map, i.e., by a plan defined by 7 = (l^d x Tp)^p. 

Let us remark that the Wasserstein distance in one dimension can be easily characterized since 
the optimal transport map, if exists, in one dimension is always a monotone non decreasing function. 
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In fact, as shown in [301 Theorem 2.18], the optimal plan in one dimension is independent of the cost 
and given in terms of the distribution functions associated to the probability measures and their 
pseudo-inverses. In fact, one can show that: 

Lemma 1.4. Given p,n (z 7'|'^(M), the optimal transport map Tp in M for d2 between them is 



2 

essentially increasing, i.e., Tp is increasing except in a p-null set. 



Proof.- Indeed, we can use the distribution function Fp(x) := p{{—oo,x)), and define its pseudo- 
inverse of Fp by the formula 

F-^{s) := sup{x G M : Fp{x) < s} , for s G [0, 1]. 

In one dimension, the optimal map for d2, and for general costs, is given by Tp = F~^ oFp satisfying 
obviously Tp{si) < Tp{s2) for all < si < S2 < 1- Thus, the optimal transport map Tp is 
nondecreasing. Since p £ 'P"'^(]R), Tp is an injective function except in a p-null set (see [21 Remark 
6.2.11]). Therefore it follows at once that Tp is increasing except in a p-null set, i.e., it is essentially 
increasing. o 

Following the seminal ideas for the porous medium equation in [28] and the linear Fokker-Planck 
equation in [23], a theory of gradient flows in the space of probability measures {V2{R'^),d2) has 
been fruitfully applied to general class of equations in the last decade [21 [9l [1] . These equations 
are continuity equations where the velocity field is given by the gradient of the variational derivative 
of an energy functional. More precisely, they are of the form 

(1.5) ^ = div (pV^) , in (0, +oo) X , 
where the free energy functional £ is given by 

(1.6) £[p]:= I U{p{x))dx+ I p{x)V {x) dx + ]- j j W {x - y) p{x) p{y) dx dy 

under the basic assumptions U : M"*" — ?> R is a density of internal energy, 1/ : M'^ — ?> M is a confinement 
potential and Vl^ : M*^ — )■ M is a symmetric interaction potential. The internal energy U should satisfy 
the following dilation condition, introduced in McCann [26] 

(1.7) A I — > X'^U{X^'^) is convex non- increasing on M"*". 

The most classical case of application, as it is for our case, is U{s) = slogs, which identifies the 
internal energy with Boltzmann's entropy. 

We can check that, at least formally, our equations of interest (jl.ip and ()1.2p are of the form ()1.5p 
in d = 1 defined by the functional (|1.6p with the choices: U = V = 0, and W{x) = log and 
U{s) = K slogs, V = 0, and W{x) = — ^ log \x\, respectively. However, the theory developed in [2] 
is not directly applicable to (jl.ip and ()1.2p for two reasons: this theory uses a convexity property 
of the functional £ that we will discuss next and the potentials V, W have to be smooth functions 
while we deal with the singular potential — ^ log |x| with an apriori unclear convexity properties. 

The needed notion of convexity for functionals on measures was introduced in [26] and named 
displacement convexity. This notion provides functionals of the form ()1.6p with a natural convex- 
ity structure allowing to show that the variational scheme introduced in p3] is convergent under 
smoothness and convexity assumptions of the confining and interaction potentials V, W and the 
convexity property of the internal energy in (jl.Tp . see [2] for precise statements. Let us define 



(1.8) Wix) 




for X 
at X = 



In our case, we will show that the interaction functional £ with U = V = and W = W given by 
p.Sp in one dimension is indeed displacement convex. The intuition behind this is that W{x) is 
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clearly convex for x > with the definition above and the optimal map between two measures is 
essentially increasing as shown in Lemma [1.4[ Therefore, when transporting measures we only "see" 
the convex part of W{x). 

On the other hand, this convexity will allow us to avoid the singularity too. In plain words, we 
will show that this interaction potential is extremely repulsive in one dimension producing that any 
initial measure is instantaneously regularized to an absolutely continuous measure for all t > 0. 
This behavior is very interesting compared to fully attractive potentials. In fact, equation ()1.5p 
has been studied in d = 1 with the displacement concave attractive potential W{x) = - log\x\, 
U{s) = K slogs and V = 0, the so-called one dimensional version of the Patlak-Keller-Segel model, 
in [7]. There, it is shown that the variational scheme in [M] converges to a weak solution of the 
equation in case the diffusion k does not go below certain critical value. Let us point out that 
nonpositive solutions to (|1.2p corresponds easily to nonnegative solutions to this ID-PKS model via 
reflection, cf. [T^ . 

Other related works for fully attractive potentials may lead to finite time blow-up, in the sense 
of finite time aggregation to Delta Dirac points, see [HIH]. For fully repulsive potentials like the one 
we consider here, we are aware about the recent work in [5] dealing with the asymptotic behavior 
of n L°°-solutions in dimensions d>2 for the Newtonian potential from a more fiuid mechanics 
perspective. 

2. Free Energy Properties 

With the purpose in mind to give a well-posedness theory for probability measure solutions to 
(jl.ip and (|1.2p . we should keep in mind that we are also interested in the asymptotic behaviour of 
the solutions. For both reasons, it is obvious that a deep preliminary study of the minimization and 
convexity properties of the free energy functionals involved has to be performed. 

In order to find self-similar solutions to ()1.2p . we will need to rescale variables, as usually done in 
nonlinear diffusion equations pT] to translate possible self-similar solutions onto stationary solutions. 
The rescaled equations can also be considered gradient fiows of certain free energy functionals which 
are uniformly 1-convex functionals in the sense of displacement convexity. These are the objectives 
of this section. 

2.1. Self-similar variables and gradient flow structure. We introduce the following self-similar 
variables 

1 



(2.1) 



y = x{l + 2t) 2 , for all t > and x G 
T = \ log(l + 2t) 



2 

Now, observe that if u{x,t) is a solution of system ()1.2p . then the function 

p{y,T) = {2t + 1)'2 u{x,t) 
with (y, r) defined by (j2.ip satisfies the equation 
(2.2) drP = dy{yp) + ndyyp - {H{p)p)y . 

Then, we can write the system (jl.2p in the new variables as 

drP = dy{yp) + Kdyyp - {H{p)p)y , Vr > and y € 
p{y,0) = uo VyGM 



(2.3) 

Equation (j2.3p has a gradient flow structure in the sense of subsection 1.2, i.e, we can rewrite it 
in the following form 



(2.4) ^ = ^ 

^ ' dt dy 



d ( 1 1 I I y 



l^/^log/j- -log|y| *p+ 2 
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where we have replaced the letter r by t again. From now on, we identify the time dependent 
probability measure p{-,t) = pt with its density with respect to Lebesgue and we use the notation 
dpt = dp{x, t) = p{x, t) dx. 

Let us begin by introducing a precise definition of a free energy functional -Ek,q on the space of 
probability measures V2iM.). We define E^^^ : 7^2(K) ^ K U {+00} as 



(2.5) 

with a,K > and where for p E V2i^) 



kU[p\ + a V[p] + W[p] for p G V^^l 
+00 otherwise 



[ p{x) log pix)dx iorpeV^^R) ^ /■ 

Jr , V[p\ ■■= / yPWaa;, 



+00 otherwise 
and y^[p] '■= 7: I W {x — y) p{x) p{y) dx dy , 

where W is given by (jl.Sp . We can now identify that ()2.3p or ()2.4p . (jl.ip and ()1.2p belong to the 
class of equations (jl.Sp with the choices W{x) = W{x), U{s) = k slogs and V = with different 
values for k and a. Thus, they are formal gradient flows of the corresponding free energy functionals 
Ef^^aip]- It can be easily checked that the functional Ef^^a is formally a Lyapunov functional for the 
equation ()2.2p . i.e. 



where 



J~-K,a[p\-= I ( ( «^logp(x) + ay - / W{x -y)p{y)dy^ p{x) dx . 



Let us remark that the functional W is also known as the logarithmic energy of p as introduced 
and deeply analyzed in |29j (see also |3H [32] ) . The next lemma shows the lower semi-continuity of 
the functionals U{p), V{p), and W(/5), and as a consequence, of the functional q,. 

Lemma 2.1. The functionals lA, V , and W are lower semi- continuous in 7^2(1^) with respect to d2- 
Moreover, the functionals lA and i?o,a with a > are weak-* lower semi- continuous in 7^2 (I^)- 

Proof. The weak-* lower semi-continuity olU is proven in [26, Lemma 3.4], which implies the 
d2 lower semicontinuity. The weak-* lower semi-continuity of V is straightforward from properties 
of weak-* sequences and it is trivially continuous for the ^2 topology. 

Before starting the proof for W, let us comment that it is essentially contained in |26[ Lemma 3.6], 
although the author deals with a more regular interaction potential W{x). The proof is inspired 
from arguments in [29l Theorem 1.3]. Here, we included it for completeness. Let us consider the 
functional 



EqAp] ■= ""'^(/o) + W(p) = / R{x, y) p{x)p{y) dx dy . 

with a > and 

[ -l-cxD for X = y 

Since the function R{x,y) — > 00 as |(a;,y)| -^00 and as |x — y| 0, it is obviously smooth except 
at the diagonal and bounded from below, then it can be approximated pointwise by an increasing 
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sequence of functions Rk{x, y) G Co°(M x M) as A: ^ oo. If />„ — )• /? weak-* as measures, then certainly 
the product measure /?„ x pn converges to p x p weak-* as measures in R x M. Now, define 

Eoa[p\-= Rk{x,y) p{x)p{y)dxdy , 

and note that E^^lpn] < Eo^a[Pn], for all n G N. Then, due to the weak-* convergence, we get 

E^Jp] = lim E^^[pn] < liminf ^o,aH, 

n— >oo ' n— ^oo 

for fixed A; G N. On the other hand, by monotone convergence Eq ^[p] £^o,a[p] as A; — )• cxd and we 
obtain EQ^a[p] < liminf„_j.oo -E'o,Q[pn]- The remaining statements follow from the continuity of V in 
the ^2 topology. o 

Remark 2.2. Let us note that the domain of the functional W consists only of absolutely continuous 
measures with respect to Lebesgue D(yV) C 'P2^{M.). This is a consequence of the definition of W 
and the fact that given a measure p with atomic or singular part in its Lebesgue decomposition, then 
p X p will charge the diagonal with positive measure. Notice that H n L°°(M) C D(W). 

2.2. Minimizing the inviscid free energy functional. The aim of this section is to make a 
summary about how to show the existence of a unique minimum among all probability measures 
in ^2(1^) to the free energy functional P[p] ■= Eq^i[p]. By Lemma |2.H we already know that P is 
weak-* lower semi-continuous, and then, in order to ensure the existence of a minimum, we only need 
to show that the functional is bounded from below. Uniqueness, compact support, characterization 
and the explicit form of the minimum of this functional were studied in relation to the logarithmic 
capacity of sets, free probability and connections to random matrices in [291 ESl ED E2]- We will 
prove some of them for the sake of the reader in the next proposition. 

Proposition 2.3. Let -0 := inf {P[p]; p e V2{M.)} . Then: 

1} -d is finite. 

ii) There is a unique p G 'P|'^(M) such that P[p] = with compact support. 

iii) Moreover, one can characterize p as the unique measure in 'P|'^(R) satisfying that 

[ fog \x - y\p{y) dy > Cp 
a.e. x G M with equality on supp{p) and with 

Cp := 2i9 - / ^p(x) dx . 

iv) Furthermore, the minimum can be explicitly computed by using the previous characterization 
and is given by the semicircle law, i.e., p is the absolutely continuous measure with respect 
to Lebesgue with density given by 



X 

y 



p{x) dx = — \/ {2 — dx . 



Proof. Part i): We first show that P[p] > for all p G 7^1^ (M) implying that 1} > —00. For this, 
observe that, for all {x,y) G M^, 



Thus, 



, _(£!±2!) ^ „ , , ,, (i^i+i.i)" ^ _£ /2V/^ ^ 

< \x — y\e 2 < {^\x\ + \y\)e 4 < supre * = - < 1- 



log ( |x — y|e ' 2° ' ] > i log(e/2) > 
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and 
(2.6) 



P[p] 



log ( \x — y\e 



> 



dp{x)dp{y) 



Uog{e/ 2)dp{x)dp{y) = Uog{e/2) > 0. 



Therefore, t9 > 0. Choose pi2 = dx\^^ 2) € 'P2'^(M), where dx denotes the Lebesgue measure. Observe 
that — log |x — y| > on (1, 2) x (1, 2). Noting that sup{|2 — y\ ,\l — y\} < 1 when y G (1, 2), then 
we obtain from Tonelli theorem that 

- / log \x - y\ d{pi2 X P12) = - / / log |x - y\dx ] dy < - / log zdz ] dy < 00 . 

JR2 J{1,2) \J{1,2) J J{1,2) \J{o,i) J 

We conclude that < +00. 

Part a): The proof follows closely |29( Theorem 1.1.3]. We start by showing that given a probabil- 
ity measure, we can always construct another compactly supported measure that lowers the energy. 
First, we observe that for a sequence with lim + |y„|) = 00, we have 

n— i-oo 

0<\xn-yn\e 2 < [\xn\ + \yn\) e 4 ^ Q, as n oo, 

and then 

lim log 

n— >oo 

Therefore, there exists sufficiently small e > 0, such that 

(2.7) - log [x - y|e-^^^ > 2(t? + 1) if (x, y) ^ x 

with := {x G M; e'^'/^ > 



+ 00. 



Next, we claim that if p G 'P(M), with supp(p) fl (M\i;e) ^ and P{p) < i? + 1, then there exists 
a /5 G ViTi^) such that P{p) < P{p)- Note this implies that there exists e > 0, such that 



(2i 



?? = inf{P(p);pG7'(S,)}. 



Thus P{p) = "!? is possible only for measures p with support in S^. 

Now, observe that ()2.7p and (|2.6p together with P{p) < + 1 implies pi^'e) > 0. This allows us 
to define 



Moreover, we have 



P[P] 



log 



\x — y\e 



1 

+ 2 



log 



(S,xS,)<= 
1 



\x — y\e 



+ 



log 



\x — y\e 



p{x)p{y)dxdy 

p{x)p{y)dxdy 
p{x)p{y)dxdy 



2{'& + l)p{x)p{y)dxdy 



2 7]R2\EeXS, 

= P[p]p(S,)2 + (^ + l)(l-p(S,)2). 

>P[p]p{^e? + P[p]{l-p{^ef), 

since P{p) <■!? + !. Hence, P[p] > P[p] and the claim follows. 
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As a consequence of (j2.8p . if p is a minimum for P, then p has compact support in S^, and thus 
p e P2(M). 

A standard argument in calculus of variations now shows that the minimum is attained in the 
set 7^2 (IK)- By definition of § and ()2.8p . there exists a minimizing sequence, i.e., {pn} ^ "^2(1^) 
with P[pn\ ^ as n ^ 00 with supp(pn) C He for all n € N. Note that each pn has support in 
the compact S^, and then we have that the minimizing sequence of measures is tight in the weak 
convergence of measures. Therefore, we can select from {pnjneN a weak* convergent subsequence 
and without loss of generality, we can assume that {pnjneN itself converges to p € 7^2 (1^) in the 
weak* topology of measures and in the 3,2 sense. Therefore, from weak—* semi-continuity of P, we 
have 

^ < P[p] < hm inf P[pn] = 

n— ^00 

and thus P[p] = i?. The absolutely continuity of the minimum is a direct consequence of the Remark 
EJ since D{W) C P|^(M). 

The uniqueness of the minimum is proven in \29[ Theorem 1.1.3]. However, in our context this 
will be clear later on from convexity properties, so we postpone this discussion. 

Parts in) and iv): The characterization of the minimum is due to the Euler-Lagrange equations of 
the variational problem with the mass constraint. However, since the minimum has compact support, 
then one obtains this variational inequality outside its support. Moreover, this characterization 
allows to find explicitly the minimum given by the absolutely continuous measure with density 
defining the semicircle law in (iv). We refer for all details to Theorem 1.1.3 and Theorem IV. 5.1] 
and |23] since it is a well-known fact in free probability and logarithmic capacity. o 

Remark 2.4. An important consequence of this result is that there exists a unique compactly sup- 
ported stationary solution of problem (|2.3|) in "^2(1^) explicitly given by the semicircle law [23j or 
the Barenblatt-Pattle profile for m = 3 of nonlinear diffusions [llj. Therefore, using the self-similar 
change of variables (|2.ip . the problem (|1.2p admits a unique, up to invariance and translations, 
global in time self-similar solution with initial data in 'P2i^)- This is already obtained and studied 
in [HIE]. 

2.3. The viscous case k > 0. In this subsection, we are concerned with the study of the functional 
Ef^^a for K > 0. Our intent is to show that functional reaches a unique minimum point on ^^'^(R). 
As we already discussed before, a suitable notion of convexity of the functional in the set of measures 
will be very important in this case. Next, we recall the definition of convexity along generalized 
geodesies of a functional E : V2{^) — > M. 

Definition 2.5. [2J A generalized geodesic connecting p to p, (with base in v and induced by 7) is a 
curve of the type gt = (vr^"^^) t ^ [0, 1], where 7 € r{i^,p,p,), vr^' 7 G To{u,p), tTj,' 7 G ro{iy,p), 
and 7r|~^^ = (1 - t)7r2 -|- tvrs. 

In particular, when dealing with absolutely continuous measures p, p.,!^ E P|'^(M) and with p = v, 
gt := [(1 — t)Id + tTp]^p is a generalized geodesic connecting p to p. In this case, we call gt the 
displacement interpolation between p and p. 

Definition 2.6. [2j A functional E : V2{^) — > (— 00, -|-oo] is \— convex along generalized geodesies 
(a.g.g. by shorten) if for every v,p,p G D{E) := {p G V20^)', E[fj] < 00} C V2(M) o,nd for every 
generalized geodesic gt connecting p to p induced by a plan 7 G T{v,p,p), the following inequality 
holds: 

E[gt] < (1 - t)E[p\ + tE[p] - h{l - t) d^^ip, p), 

where 

dlip,p)-= / \x3 - X2f d'y{xi,X2,X3) > dl{p,p) . 
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If gt is the displacement interpolation and X = 0, we say that the functional E is displacement convex 
as originally introduced in [26j. 

We readily apply these notions of convexity to our functional E^^a- 

Proposition 2.7. Let K,a >0. The functional E^^a defined by (j2.5p is a— convex along generalized 
geodesies. 

Proof. Following the notation in (|2.5|) . we can reduce ourselves to show that the functional 
W is convex along generalized geodesies in D{E^^a) C Vif'iM.). In fact, it is well-known that lA is 
convex (A = 0) along generalized geodesies, and that the functional V is 1— convex along generalized 
geodesies in D{E^^a) C 7'f'=(M), see [2l[26j for details. 

Let p, /U,i/ G D{E^^a) C ^'^'^(M) and gt be a generalized geodesic connecting p to fj, with base 
point v. As we are working on the real line M, we can express the generalized geodesies as 

gt = ((1 - t)TP + tTlf)#iy, 

where Ty and Tlf are the optimal transport between v and p, and u and p respectively, with the 
properties in Lemma ll.4[ 

Let us observe that by definition of gt and its absolute continuity with respect to Lebesgue, we 
get 

W[gt] = / -\og(\x -y\)d{gt x gt) 

= / -\og{\{l -t){T{:{x)-TP{y))+t{Tlf{x)-Tif{y))\)d{uxu) 

<{l-t) I -\og{\{TP{x)-TP{y))\)d{uxu)+t [ - log{\{T>f {x) - Tjf {y))]) d{u x u) 
= {l-t)W[p]+tW[p], 

where the last step follows from convexity of the function — log |a;|. This last step is fully rigorous 
provided the following claim (C) holds: there exists a i^-null set A such that 

(1 - t)iTPix) - TP{y)) + tiTlfix) - Tlf{y)) + 0, 

for all t G [0, 1], x,y and x ^ y. In other words, the interpolation map reaches the logarithmic 
singularity only if x = y or in a x z^-null set. 

In order to prove this claim, we remind that the optimal transport on the real line between two 
measures in 1^2^ {M.) is essentially increasing, see Lemma ll.4[ Now, let A be a i/-null set such that 
Tj^ and are increasing in A^. If x^y ^ A'^ and x ^ y then let us show that 

(1 - t){TP{x) - TP{y)) + t{Tlf{x) - Tlf{y)) / 0, Vt G [0, 1]. 

To prove this, suppose that 3 to £ (0, 1], x, y G yl^ and x ^ y such that 

(2.9) (1 - to){TP{x) - TP{y)) + t,{T>f{x) - Tjfiy)) = , 

then we deduce 

{T>fix)-Tlf{y)) to-1 

{TS{x)-TS{y)) to - ' 
which provides a contradiction, because the optimal transport maps Tlf and are increasing in 
A^. For the case to = in (|2.9p . we have TS{x) — Ty{y) = which yields x = y or x, y G ^4, because 
of essentially injectivity of Ty . This finally shows the claim (C), and thus that W[5't] is a convex 
function in t G [0, 1] for all generalized geodesies corresponding to absolutely continuous measures, 
which gives by definition the convexity of W along generalized geodesies in D{Ei^^a)- o 

Proposition 2.8. Let k > 0, a > 0, and -d^^a ■= inf {E'k.qIp); P G 'P2(^)}- Then: 
i) '&K,a is finite. 
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ii) There is a unique PK,a ^ 1-^2^ {^) such that EK.^a\pn,a a ■ 

Proof. This results is proven for k = in Proposition 12.31 except the uniqueness part. Let us 
assume from now that k > 0. 

Part i): Recalhng the definition of the functional in terms oiU, V, and W we split 

a 

It is straightforward to use Jensen's inequality to show that 

^[P] + f V[p] = log (plv^) e-'^^'/'dx > |(log5 - log(2vr)), Vp G VI 

with 5 = a 12k. Proceeding analogously to the proof of Part i) of Proposition 12. 3| we obtain 

_ (\x\ + \y\f _ ri / 4 \ 

< |x - y|e 4 < (|x| + \y\)e " s < supre s = — , 

r>o \aej 

and 

(2.10) ^o,a/2[p] = - log ( - J dp(x)dp{y) > -- log i—j . 

Therefore, the functional -E^jO is bounded from below, and so > — oo. Since the domain 

D{Ef^^a) 7^ 0) there is /o € 7^2(1^) such that Ei^^a{p) < oo, and therefore the infimum is finite. 

Part ii): From Lemma 12.11 is weak-* semi-continuous in 7-'f'^(M), which shows that the 

infimum is achieved at some point p € 7^|'^(M). In fact, it is easy to check based on the same 
arguments for Part i) and Proposition 12.31 that any minimizing sequence is weakly compact in 
-L^(M), see similar arguments in [7J, since k > 0. 

The uniqueness claim for k ^ follows from the strict displacement convexity of Ei^ q^ . Indeed, let 
pi and p2 be two different minimum in "^(M) to -E^.a and consider pi the displacement interpolation 
between pi and p2 sd, t = 1/2. By the a-displacement convexity of Ef^^a, we have 

which provides a contradiction. Therefore, there exists a unique minimum of * 
3. WeLL-POSEDNESS, ASYMPTOTIC BEHAVIOR AND INVISCID LIMIT 

As pointed out in the introduction we will obtain solutions for ()2.3p as the limit of a Euler 
approximation scheme in probability space V2{^)- More precisely, consider a step time r > and a 
initial data po ^ 'P2{^)- We define, for a fixed p, the functional X(r, p, .) : (IK) — > (-co, oo] as 

I{t,p;p) := ^dl{p, p) + E[p]. 

Formally, we define the following recursive sequence (p")„gN: 

(3.1) p° := po 

(3.2) p^= mm I {t,p^~';p), n e N, 

pGP2{K) 

which can be seen as the discrete approximate Euler solution to the gradient flux equation 

^ = -VE[p], t>0 

in the metric space (V2(^), ^2)- More precisely, one calls a discrete solution, the curve pi obtained 
as the time interpolation of the discrete scheme (|3.ip - (|3.2p connecting every pair (p"~^,p!^) with a 
velocity constant geodesic in t € [(n — l)r, nr), see |2J. 
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3.1. Gradient flows. Below we remember the definition of a gradient flow solution. 

Definition 3.1. We say that a map pt € ACi^^{{0,oo);V2{^)) is a solution of the gradient flow 
equation 

-vt G dE{pt), t > 0, 
if —Vt G Tanp^V2i^) belongs to the suhdifferential of E at pt, a.e. t > 0. 

It is known that pt being a gradient flow in V2 (IR) is equivalent to the existence of a velocity vector 
field —Vt G Tanp^V20^)^\ dE{pt) a.e. t > 0, such that ||t't||j^2 € L^^^(0, 00) and the continuity 
equation holds in the distribution sense: 

(3.3) + V .{p^vt) = Q in Mx(0,oo). 

The next theorem ensures the existence of a gradient flow solution for the free energy functional 
Ef^^a as in 



Theorem 3.2. Let k, a > 0, /?o G 7^2(1^) o-nd the functional E^^^. The following assertions hold: 

(1) (Existence and Uniqueness) The discrete solution pi converges locally uniformly to a 
locally Lipschitz curve pt := >S't[/Oo] in 7^2(1^) which is the unique gradient flow of E^^a with 
limt_^o+ Pt = Po- Moreover, the curve lies in P|''^(]R), for all t > 0. 

(2) (Contractive semigroup) The map t 1-^ 5't[po] for all a > is a a-contracting semigroup 
on i-e- 

d2{St[pQ\, St[pQ]) < e~"*(i2(po,/Wo) for all po,Mo G ^2(1^) • 

(3) (Asymptotic behavior) Let q > and let us denote by p^^a the unique minimum of E^^^^. 
Then for all < to < t < 00, we have 

ci2(pt,P«,J<e-"(*-*«)d2(pi„,p,,J 

and 

E^M - E.AV.A < e-'"(*-*«)(i^.,a[pto] - • 

(4) (Free Energy identity) The solution pt := St[po\ is a curve of maximal slope and it satisfies 
the identity: 



En.,a[Pt\= EnAPs\+ / / \Vr{x)\^dxdT 
J s JR 



for alio < s < t, where vt G Lf^ci^, 00; L^^(]R)) is the associated velocity field satisfying 
in Detinition \3.1\ . 



Proof. First notice that T'2nL^nL~ C D{E^^a) C 7'f''(K) by Remark [2:21 and D{E^^c,) = P2(K). 
Let us start with the case a > 0. Collecting the results obtained through previous sections, we 
have that the functional E'k,^ : 7^2(1^) (—00, +00] is a proper, l.s.c, coercive functional and 
a— convex along generalized geodesies. Moreover, X(^, r; p) admits at least a minimum point /i,-, for 
all T G (0, T^,) and jj, G ^2(1^)- The minimum p^ G ^2*^(1^) for k > 0, because minimizing sequences 
are weakly compact in L^(M) as in Proposition 12. 81 and by the Remark l2.2l if k = 0. Therefore, all the 
statements result directly from the general theory of gradient flows developed in [21 Theorem 11.2.1]. 
In case a = 0, we deal with plain convex functionals along generalized geodesies and the same results 
apply. However, we need to be careful with the coercivity and the existence of minimizers for the 
one-step variational scheme since we lack of a direct confinement. This is easily provided by the 
following observation using the triangular inequality and 

(3.4) l{T,p-p) := l-4{p,p) + E^4p] > l-dl{p,5o)-^dli6o,p) + E^^o[p] = -^Vip] + E^^2/r[p] 
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for all p (z 7^2 (IR) and all r > 0. Therefore, this implies the boundedness from below, the existence 
of minimizers of the one-step variational scheme, and the coercivity in the case of a = 0. Again the 
results of p] Theorem 11.2.1] apply directly. o 

Remark 3.3. As a consequence of previous theorem, we have shown the global-in-time well-posedness 
for the Cauchy problem for general measures in V2 (M) as initial data for equations (jl.ip and ()1.2p 
and their self-similar counterparts ()2.3p . Moreover, we have shown the convergence towards self- 
similarity in the sense expressed in the third part of Theorem 13. 2L Note that the gradient flows 
obtained in Theorem \3.2\ for the functionals Ei^^ and E'^^i are equivalent through the change of vari- 
ables ()2.ip . Finally, note that the evolution is defined in a unique way for any initial data in 7^2 (I^)- 
However, the evolution flow regularizes instantaneously since it belongs to 'P|^(M) for allt > 0. This 
is the precise mathematical statement showing that the repulsive logarithmic interaction potential in 
one dimension is "very repulsive" . 

Remark 3.4. (Power-law potentials) Consider the power-law interaction potential W{x) = \x\~^ , 
for < /3 < 1, and its natural extension 



(3.5) Wsix) 




for x ^ 
at X = 



Let the free energy functional Ep^i^^^ '■ "^2(1^) — )• M U {+00} be defined as in ()2.5p with Wp instead of 
W. Noting that \x\~^ is locally integral when < /3 < 1, similar arguments as in Propositions 12.31 
and \2.S\ give that the infimum 'dp^^^a is finite and unique. Moreover, 7^2 D D (R) C D{Ep^^^a) C 
■P2"^(M) and = 'P2{R)- Since the function \x\^^ is convex, we have that -E^^^.a is a— convex 

along generalized geodesies. Thus, the results on gradient flows given by Theorem \3.2\ also hold true 
for ()3.5p and we can take initial measure pQ € 7^2 (I^)- 

For range f3 > 1, we have D{W) = {0} and D(Ei3^f^^a) = {0} and then the theory trivializes. In 
particular, if p € D{W) were a continuous function and if = {x G M : p{x) > -} is a positive 
measure set, for some n € N, then 

If 1 

00 = — / adxdy < W(/9) < 00, 

n'^ JkxK \x - yf 

which gives a contradiction. Therefore, the only continuous function in the domain is p = 0. 

Theorem 3.5. (Inviscid limit) Let us consider the functionals and E^^a with a > 0, corre- 
sponding to viscosity k > and k = respectively, and assume that pQ G D{E^^^a) with eo > 0. // 
Pnii), p{t) are the corresponding gradient flow solutions in 7^2 (IK) with initial data po, then 

p^{t) ^ p{t) m T2{R) 

locally uniformly in [0, 00), as k — )■ 0+. 

Proof. In view of the stability property of [21 Theorem 12.2.1], we need only to verify (in a 
neighborhood of k = 0) the equicoercivity of the family of functionals {-E're,o}K>o and the uniform 
boundedness at pQ. More precisely, we need to show 

sup Ef, a[po]<'X> and inf -^djip, p) -\- Ef, a[p] > -00, 

Kg{o,£o) «e(o,eo), P6P2{K) 2r 

for some r > and p £ 7^2 (1^) and eo > 0. Firstly, observe that po € D{E^g^a) implies U[po] < 00, 
V[po] < 00 and W[po] < 00. It follows from ([23]) that po G D{Ef^^a) for ah k G (0, 00). Also, 

sup -EK,a[po] < max{0, eoZ^[/Oo]} + aV[pQ] + W[pq] < 00 . 

KG{0,eo) 
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In order to conclude the proof, it remains verify the equicoercivity. By using ()3.4p . we observe 

^4{^i, p) + E^^^[p] > ^dlip, p) + E^^o[p] > -^V[p] + E^^y2r[p] 

for all Q > 0. Let us split the functional i^'K.i/r &s 

E^,l/2r[p] = f^[p] + ^V[p] + W[p] = i^0,l/4r W + td^[p] + ^V[p] . 

Let us now remark that kU + is bounded from below. Note that it is the relative logarith- 
mic entropy functional leading to the classical linear Fokker-Planck equation whose minimum is a 
Gaussian M{x) determined by 

M(x) = (4vr-) exp^- — 

Therefore, we get 

(3.6) KU[p] + ^V[p] > kU[M] + ^V[M] = -\\og (47r-) > - J log (4^^ 

for K G (0,eo). Due to Proposition 12.81 and ()2.10p . -E'o,i/4r is also bounded from below by '!9o.i/4r- 
Using ()3.6p with a = 27, we conclude 

. ^ [^dl{p,p) + E^^o.[p\ \ > --V[p] - J log (Svreor) + i?o,i/4r • 



3.2. Solutions in sense of distributions. An important point is to know whether the gradient 
flow solutions are solutions in the sense of distributions. Firstly, let us define the notion of weak 
solution which we deal with. We say that a measure pt is a weak solution to equation (j2.3p . with 
initial condition pQ, if for all <p € Cq- 



(3.7) 4 / ^(.^)dPt = 1^ f ^"{x)dpt - [ ^'{x)xdpt + \ ( V^'(^) "^'^y^ dipt X Pt) 



dt Jk 2 J^2 X — y 

in the distributional sense in (0, 00) with pt po weakly-* as measures. 

In order to obtain a connection between gradient flows and weak solutions, we need to describe, 
see [2], the minimal selection of the subdifferential of -E = Ef^^i, that is the set d° E^^i{pt). For 
that matter we need to consider the following functional: for each fixed p G 7^2 (IS) we define 
L„ : CoHK) ^ M as 



p ■ ^0 
(3.8) 

Lp(^)= lim / —^^(x)d(pxp)= lim - f '^C^) " '^(^) ^(p x p) < oo,Vy g C). 
S^Q+ J\x-y\>5 X - y ^ J\x-y\>5 x-y 

It is straightforward to check that 

(3-9) \Lp{v)\ < \ ||<^'|L Ap ^p) = \ Ik'lloo ' 

and therefore Lp G (Cg(M))*. So, there exists p G A^(M) and a constant cq such that (see [121 p. 225]) 

Lp{ip)= [ if'dp + coifiO) V^gCqHm). 
Jr 

From (j3.9p . we can also see that cq = and we obtain the following representation to Lp-. 

(3.10) Lp{ip) = [ if' dp, G C^iR). 

Jr 



MASS-TRANSPORTATION APPROACH TO A ID MODEL WITH NONLOCAL VELOCITY 15 

Lemma 3.6. Let k > and n as mentioned above. If a measure p € D{Ei^^a) C V20^) belongs to 
D{\dEK^a\) then we have {kp + p) € and 

(3.11) puj = dx{Kp + p) -\- apx for some uj G Lp(]R). 

In this case the vector u defined by ()3.1ip is the minimal selection in dE^^a[p], i-s. uj = d°Ei^^a[p]- 

Proof. For each compactly supported smooth test function ip G C^(M), let us consider the map 
tpe := Id + eip. It is easy check that € D^E^^a), when p e D^E^^a) and for e > small enough. 

Since p € D{\dEi^^a\) and from the definition of metric slope see [2], we have 



\Ai{ip) + A2iip) + A3{^)\ :-- 



^ U[A*p]-m ^ ^ V[A*p]-V[p] ^ W[A*P] - W[p] 



< \dE^a\[p\ hm ^ < oo. 

e— !>0 e 

The terms Ai and A2 can be exactly treated as in [21 Chapter 11] and one obtains 

Ai{^) = —K / If' dp and A2{(p) = a xipdp. 



Now, we deal with the term A3. Notice that the map 

X y)= i -iogl(^-?/ + g(v^(^) -"pjyM -(-iog|^-y|) 

2 s 

is nondecreasing in e > 0, for fixed x,y ^ 0. As Ai{{p) and A2{ip) are finite, then ^3(9?) is also finite. 
By the monotone convergence theorem, we have 

Asiv) = hm / Qie, x, y)d{p x p) = [ ^^^^^^^dip x p) 
= —Lp{ip) = — Lf' dp < 00. 



Notice that the second integral above is not a singular integral because ip £ Cg°(M). Observing that 



e^O e 



we get 



{Ai + A2 + As){(p) = / -Kip'{x)dp + a / XLp{x)dp - Lp{tp) 
Jm. Jr 



Changing (p by —(p, we finally obtain 



{-Kp' + axip)dp - Lp{p) < \dEi^^a\[p]\\^\\Lj{ 
So, there exists uj G L'^i^) with ||<^||^2(ig) < such that 

ujpdp = {Ai + A2 + A3){ip) 
(3.12) ~ ~ { <p' dp + Lp{ip) \ + a xp dp, G C( 
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Thus u) G dE^^a[p\ and uj is the minimal selection in dE^^a[p\, i.e. a; = d°Et^^a[p\. Finally, let us 
characterize uj. Since p € ^2(1^) implies that i^lf] = f^xipdp is bounded in Lp(M) (with norm at 
most J^x"^ dp), we get 



\< dx{Kp + p),Lp >\ 



ip' d{np + ji) 

< (^\dE^,M + a l^x^dp^ . 
Therefore, dxi^p + /i) £ AA{M.), i.e. (k/? + /i) € BV(R). Integration by parts holds: 



ip{x)d{dx{Kp + p)) 



m\Li{ 



which implies d^inp + p) e L'j,{R)nM{R) and then {np + p) e W^^ 



1,1 



Finally, coming back to 



(|3.12p . we obtain the following expression for uj, the element of minimal norm in the subdifferential 
of Ef^^a'- P^ = dx^Kp + p) + ax p. o 

The next theorem gives a connection between gradient flows and the notion of weak solution ()3.7p . 

Theorem 3.7. (Distributional solution) Let pt correspond to pt through \3.^) . For every po G 
7-'2(IK) and every k, q > 0, the gradient flow pt in 7^2 (1^) of the functional E^^a is a distributional 
solution of the equation 



(3.13) 



dpt 
dt 



d_ 

dx 



d_ 

dx 



Pt 



dxiK-Pt + Pt) 



Pt 



+ ax 



satisfying p{t) po as t 0+ , pt e Ll^^ii^, +00); W^^^ (M)), and 



d_ / Oxjupt + Pt) 
dx V Pt 



+ ax 



e^ioc(o,+oo) 



Proof. Since ujt = d° E^^^ipt) and pt is the gradient flow of E^^^^ it follows from Lemma [3.61 that 
(j3.13p is satisfied by pt with the additional conditions found in the statement of the theorem. 
Now, observe that pt satisfies (|3.13p is equivalent to pt satisfies 

dec c 

ip'{x)d{dx{Kpt + Pt)) - a / ip'{x)xdpt 



dt 



ip{x)dpt 



ip"{x)d{Kpt + Pt) - a I ip'{x)xdpt 



K I ip"{x)dpt - a I (p'{x)xdpt + / ip"{x)dpt 

Lp'{x) - cp'iy) 



= K (p"{x)dpt - a <p'{x)xdpt + - 
Jr Jr ^ 

in the sense of distributions on t € (0, oo) and for all 99 E Cg* 



x-y 



d{pt X Pt), 



Remark 3.8. Although in the general theory in [2\ the previous result is a characterization of 
gradient flow solutions, we do not know how to get the converse in the characterization of the 
element of minimal norm in Lemma [3.61 since we do not how to show that p is absolutely continuous 
with respect to p. This implies that we do not know how to show that distributional solutions with 
the properties written in Theorem 13.71 are gradient flow solutions. 

Remark 3.9. (Power-law potential) In the case of potential W{x) = \x\~^ with < /? < 1, we also 
have that the associated gradient flows (see Remark \2>A\j is a solution in sense of distributions as in 
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(|3.7p or ()3.13p . Instead of ()3.8p . m i/izs izme i/ie operator Lp{ip) is given by 
Lpiif) = lim / _/3 J^^|) ^(x)d(p X p) 

= lim i / -S ' /'"'-^fa> rf(p X p) < CO, 

/or (/3 G Co(M). r/ius, since p € D{Efs,K,a), 

^ ./r2 |x — yf 



for all if G Cq(M). Therefore Lp G (Co(M))* and, similarly to ()3.10p . Lp{ip) = J^ip'djjL, for ip G 
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